Classical and relativistic node precessional effects in 
WASP-33b and perspectives for detecting them 
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Abstract WASP-33 is a fast rotating, main sequence 
star which hosts a hot Jupiter moving along a ret- 
rograde and almost polar orbit with semi-major axis 
a = 0.02 au and eccentricity provisionally set to e = 0. 
The quadrupole mass moment and the proper angu- 
lar momentum Si, of the star are 1900 and 400 times, 
respectively, larger than those of the Sun. Thus, huge 
classical and general relativistic non-Keplerian orbital 
effects should take place in such a system. In particular, 
the large inclination of the orbit of WASP-33b to the 
star's equator allows to consider the node precession 17 
and the related time variation dtd/dt of the transit du- 
ration td- The WASP-33b node rate due to is 9 x 10^ 
times larger than the same effect for Mercury induced 
by the Sun's oblateness, while the general relativistic 
gravitomagnetic node precession is 3 x 10^ times larger 
than the Lense-Thirring effect for Mercury due to the 
Sun's rotation. We also consider the effect of the cen- 
trifugal oblateness of the planet itself and of a putative 
distant third body X. The magnitudes of the induced 
time change in the transit duration are of the order 
of 3 X 10-6,2 X 10"'', 8 X 10"^ for J^, the planet's ro- 
tational oblateness and general relativity, respectively. 
A yet undiscovered planet X with the mass of Jupiter 
orbiting at more than 1 au would induce a transit du- 
ration variation of less than 4 x lO"^. A conservative 
evaluation of the accuracy in measuring dtd/dt over 10 
yr points towards « 10~®. The analysis presented here 
will be applicable also to other exoplanets with similar 
features if and when they will ne discovered. 
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1 Introduction 



HD1 5082 (WASP-33) is 



an early type, A5mA8F4 
located at 116 pc from us 
( Perrvman et al.l |1997() , which has the peculiarity of 



star ( Grenier et al. 19991) 



hosting a close hot Jupiter moving along a 1.22 d ret- 
rograde orbi t highly inclined with respe ct to the stel- 



lar equator ( Collier Cameron et al. 2010f ). Other plan 



ets with the same char acteristics have been recentl y 
discovered: WASP-2b dCollier Cameron et al.l boOTi) . 

2008t 



WASP-5b ( Anderson et al 



Gillon et al 



200SL 



WAS P-8b dOueloz et al.ll2010l). WA SP-15b (IWest et al 



20091). WASP-17b (lAnderson et al.ll2010l) and HAT-P 



7b iNarita et al.ll2009t IWinn et al 



2009f ). Such orbital 



features are unexpected since the stars and the plan- 
ets usually originate from the same disk. A possible 
explanation may be the presence of yet undiscovered 
planets which have acted in the past in such a way 
to induce the misalignment of the systems' angular 
momenta ( Anderson et al.lboiot Collier Cameron et al 
2010l ). Evidence for a distant third b ody in the HAT-P - 



7 system has been, actually, found ( Winn et al. 20091) 



also WASP-5 may host a third body ([Anderson et al. 
20081) . 

WASP-33 rotates rapidly, almost 1 — 2 orders of mag- 
nitude faster than the other stars hosting planets mov- 
ing along highly inclined orbits, so that it is expected 
that its shape presents relatively large deviations from 
sphericity. Indeed, as we will show, its quadrupolar 
mass moment should be 1900 times larger than 
that of the Sun, while a reasonable evaluation of the 
magnitude of its proper angular momentum 5*^ points 
toward a value about 400 times larger than the so- 
lar one. Such figures, in addition to the small orbit 
of i ts planet, having a semi- maj or axis a = 0.02555 



au (jCollier Cameron et al.l I2010I ). suggest that huge 



non-Keplerian orbital effects should take place in the 
WASP-33 system. They are larger than those occurring 
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in the other similar planetary systems because of their 
wider orbitfQ (a = 0.032 - 0.079 au). Such effects could 
be fruitfully used to determine or constrain some key 
stellar and/or planetary parameters. To this aim, let us 
note that if, on the one hand, the periastron to cannot 
be used because of the provisio nally assumed circular 



ity o f the orbit of WASP-33b (|Collier Cameron et al 
2010l) . on the other hand its high inclination to the 



stellar equator makes feasible, at least in principle, the 
use of the node 17 which may translate into observable 
effects of the ph otometry of the system like, e.g., the 
transit duration (|Miralda-Escudell2002l) . 

Among the dynamical effects considered, we will ex- 
amine also the im pact of the gener al relativistic grav- 
itomagnetic field ([Mashhoonl 120071) , induced by bod- 
ies' rotation, which in the test partic le limit yields the 
Lense-Thirring effect ([Lense fc Thirring:,1918.) . Indeed, 
it is the only general relativistic effect affecting the 
node of an orbiting body, and its expected precession 
for WASP-33b is 3 x 10^ times larger than the Lense- 
Thirring node rate of Mercury in the field of the Sun. 
Moreover, as we have noted before, it is likely that 
WASP-33b is not the only planet hosted by HD15082; 
if other distant bodies will be actually discovered, it 
will be crucial reconstructing the dynamical history of 
the entire planetary system. In doing that, even small 
accelerations like the gravitomagnetic one, although 

4 X 10~^ times smaller than the dominant Newtonian 
monopole, may play a non negligible role. Such a view 
is supported by the fact that the inclusion of general 
relativitjH in the latest long-term numerical integra- 
tions of the equations of motion of the solar system 
( Laskar fc Gastineaul 2009[ l has drastically changed the 
probability of the occurrence of a secular resonance be- 
tween Mercury and Jupiter, which would induce a fatal 
increase of the eccentricity of the orbit of Mercury, from 
some tens percent to 1%; the general relativistic accel- 
erations of Mercury and Jupiter are just 7 x 10~^ and 

5 X 10~^ of the corresponding Newtonian monopoles, re- 
spectively. Last but not least, a gravitomagnetic effect 
may alias the signature of a yet undiscovered further 
planet X for certain values of its mass and distance. 

The influence of some rotation-independent, grav- 
itoelectric effects of general relativi ty in some ex- 
trasol ar planets has been t reated in Miralda-Escud^ 

" lorig (,2006, ): 



( 20021) : lAdams h Lauehlinl (l2006allbl 



^Actually, WASP-5 harbors a p lanet with semi-major axis a = 
0.027 au jAnderson et al]|2008li . but the rotation of the star is 
21 times slower than that of WASP-33. 

^Actually, in this case the general relativistic terms considered 
are those which depend only on the masses of the bodies, dubbed 
gravitoelectric l|MashhoonH2007t) . not the gravitomagnetic ones. 



Hcvl fc Gladman (.20071): [Jordan &: Bakosl(l2008l) : [Pal fc Kocsi^ 
(|2008l) : iRagQzzine fc Wolj (l2009l ): iLiI (I2OIOI) . In par- 
ticular, the perspectives in measuring the gravitoelec- 
tric correc tion to the t hird Kepler law have been con - 
sidered bv llorid (I2OO6I ) and iRagozzine fc Wold (l2009l). 

Miralda-Escudi(l2002l) . lHevl fc GladmanI (120071) .Tjordan fc Bakosj 



2008 ). iPal fc Kocsid (120081 ) and IRagozzine fc WoU 



20091) dealt with the possibility of detect ing th e gravi 



toelectric periastron precessions, while iLi (I2OIOI) looked 
at the gravitoelectric secular change of the mean 
anomaly connect ed to the variation of the periastron 



time transit t 



Adams fc LaughlinI |2006allbl l3) stud 



ied the impact of the general relativistic gravitoelec- 
tric terms on the long-term, secular interactions among 
multiple planetary systems. 

The paper is organized as follows. In Section [5] we 
set the scene by evaluating some physical parameters 
of the star and the planet which enter the dynamical 
effects analyzed subsequently. In Section[3]we work out 
the node precessions due to the oblateness of the star 
and the planet, a putative third body X, and general 
relativity. Section 2] summarizes our conclusions. 



2 Setting the scene: the relevant parameters of 
the system 

In Table [T] we quote the relevant parameters of the host 
star. A careful discussion is required about the value 
of the star's angular speed S^,, evaluated from 



Ri, sin 



(1) 



quoted in T a.ble[Il In eq. dH Uj, = t'^.sin /^ = 86.48 ± 
0.06 km s~^ ( Collier Cameron et aLllioiol ). and is the 



angle between the stellar spin axis and the line of sight, 
not to be confused with the ang le i ~ 87.67 ± 1.81 
deg ( Collier Cameron et al. 2010f ) between the orbital 
angular momentum and the line of sight. Both such 
angles are defined in such a way that = i = when 
the angular momenta and the line of sight, oriented 
towards the observer, are parallel, and = i = tt 
when they are antiparallel. Actually, is unknown; 
a direct determination of it may be possible by mea- 
suring and int erpreting asteroseismological oscillations 
(jGizon fc Sola nki 200 3i) , photometric mod ulations pro- 



duced by starsDo tj^ (iHenrv fc Winn 



20081) and polar- 
ized light curves (jFluri fc Berdvuginall2010l ). Anyway, 



•^The Kepler mission would be able to accomplish such tasks, but, 
unfortunately, it cannot see WASP-33. Anyway, time has been 
awarded on the MOST satellite to analyze the transits, search for 
eclipses and characterize the pulsation spectrum of the host star. 
(Collier Cameron A., private communication, Juno 2010). 
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Table 1 Relevant stellar parameters (jCollier Cameron et al.ll2010l ). We used M© = 1.98895 x 10^° kg, Rq = 695.99 x 10** 
m. In deriving the star's angular speed H* = sin 7* ), with = Vj, sin Ij, , see the discussion i n the text. The value 

displayed for the tidal parameter q* is an upper bound corresponding to (|Collier Cameron et al.ll201(|) m'^„^ = 4.1 mjup for 
the p lanet's mass rrip. The value of the relative semi-m ajor axis used in q* is a = 0.02555±0.00017 au (jCoUier Cameron et all 
I2OIOI ). The figure for J| has been obtained by using (Claretlll995l ) fcj = 0.03 for the star's Love number. The moment of 
inertia factor is defined as a* = C*/M*7?*q, so that the star's proper angular momentum is — a*M<,7?JqH*. The value 
used for the flattening f^, is the average of /^j^ and /max- 



Parameter (units) 


Numerical value 


(Mo) 


1.495 ±0.031 




1.444 ±0.034 


Ui, (km s""'") 


86.48 ±0.06 


2. (s-i) 


8.6 X 10-5 




0.0379461 


it 


-1.42772 X 10-"^ 




3.80175 X lO-'^ 


h 


0.0302587 


a* 


0.277011 


(kg m^ s-i) 


7.16483 X 10"^^ 



some reasonable, a-priori constraints can be posed on 
li,. First of all, the natural condition 



so that 



sin/* > 

c 



(2) 



must be satisfied in order to keep the linear rotation 
speed V* of the star smaller than the speed of light c. 
It implies a preliminary range of admissible values 



0.016 deg, 
179.983 deg. 



(3) 



In fact, tighter constraints on /* come from the con- 
dition that the centrifugal acceleration at the equator 
must be smaller than the gravitational acceleration at 
the star's surface, i.e. it must be 



It implies 
i:^,^ - 11.22 deg, 

/^,, = 168.77 deg. 



(4) 



(5) 



By numerically investigating the star's angular speed 
S* as a function of /* according to eq. ([T]) within the 
range of eq. ([5]), it turns out that 



8.6 X 10-5 s-i, 
4.5 X 10-4 s-\ 



(6) 



5.2. 



(7) 



The minimum occurs for 7^, = 90 deg. Anyway, it re- 
sults that Sj, remains below 1.5 X 10-"* s-i in about 70% 
of the allowed range for 7*; thus, we conclude that SJ^jj^ 
is somewhat representative of the most likely values for 
S. In fact, such an argument can be made even more 
plausible and stringent by numerically investigating the 
ratio of the stellar centrifugal acceleration to the gravi- 
tational one as a function of the allowed values of 7^,. It 
turns out that the tails in which S* > 1.5x10"'' s~'^ cor- 
respond to values of the star's centrifugal acceleration 
larger than 20% of the gravitational pull. Such figures 
are clearly highly unrealistic in view of the unobserved 
associated surface phenomena which, instead, would 
take place. Just for a comparison, Af^^/A'^^^^ sa 10"^ 
for the Sun. It results that for 7* = 7* the star's cen- 
trifugal acceleration reduces to about 4% of the gravi- 
tational one, which is its minimum; the same substan- 
tially holds for about 60 deg ^7* < 120 deg in which 
S*/SJ^jjj « 1.1. The previous arguments motivate our 
choice of 'E.^-^ for the stellar angular speed in Table [T] 
and in the rest of the paper. 

Concerning the first even zonal harmoni c J2 , it has 
been computed as iRagozzine fc Wolj|2009l) 



^2 



(8) 
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whercEl (|Ragozzine fc Wolj|2009l ) 



9r = 



(9) 



By comparison, the Sun's proper angular mo mentum, 
obtained fr om asteroseismology, amounts to ( Piiperd 
19981120031 ) 



50 = 1.9 X 10^1 kg m^ s 



are related to the rotat ional and tidal potentials of the 
star, and |Claretl fl995l) fc2 ~ 0.03 for mam sequence 
Sun-like stars. From Table [T] it results that the contri- 
bution to J2 of the tidal distortion raised by the planet 



(17) 



i.e., it is about 377 times smaller than eq. (|16|) . 

In Table [2] we quote some relevant parameters of 
the planet. We have assumed that the planet's an- 



of mass mp on the star is negligible. It is interesting gular speed equals its orbital frequency n; for the 



to note that the value of Table [T] for is about 1900 
times larger than t hat of the Sun, whic h is of the order 



of Jf = 2x lO"'^ (jPireaux et al.ll2007h . 



In regards as the flattening, defined as 



planet's Love number we assumed the moderate valuqj 
fc^ = 0.3 (|Ragozzine fc Woljliooj ). Note that, contrary 
to the star case, the tidal parameter 



^cq ^pol 



(10) 



9t 



(18) 



cq 



is about 3 times larger than the rotational parameter 



it is related to the adimcnsional moment of inertia pa- 
rame ter a^, by the Darwin- Radau relation ( Murrav fc DermQt)t^ _ ^p'^p 
200d) ' 



The obvious requirement that > yields 

/* > /min = = 0.0130849, (12) 

while the the condition that < 2/5 = 0.4, valid only 
for a homogeneous spherical body, implied 



(19) 



i.e. the tidal bulge raised on the planet by the star is 
larger than the centrifugal one due to its rotation. As 
a consequence, the tidal contribution to dominates 
that due to the rotation; indeed, 



fgP = 8.38747 xlO^^ 
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qf = 1.25782 xlO"^ (20) 



Concerning the planet's proper angular momentum, re- 
peating the same reasonings as for the star yields 



A < /:.ax = = 0.047326. 



(13) Sp 1.47243 X 10^^ kg m^ s'^ 



As a reasonable figure for the star's flattening we take 
the average 



f*. + f* 

J min ' J mi 



= 0.0302587, 



which yields 



0.277011, 



(14) 



(15) 



of the same orde r of magnitude of the evaluations by 
Ford et al.l (|l999l ) for main-sequence stars. As a conse- 
quence, a realistic evaluation for the star's proper an- 
gular momentum is 



= 7.16483 X 10''^ kg m^ s 



(16) 



(21) 



which is about 48000 times smaller than S*^. 



3 The orbital effects 

Since the orbit of the pla net is assumed circular 
( Collier Cameron et al.ll2010l ). it is not possible to con- 
sider, not even in principle, either the periastron u or 
the time of passage at periastron tp. Thus, we will not 
investigate the precessions induced in w and the mean 
anomaljQ A^. Thus, let us consider the longitude of the 
ascending node whose secular change is related to the 
precession of the orbital plane. Such an effect may im- 
pact some directly observable quantities in photometry 
of transiting extrasolar planets like the transit time 



^Actually, in the relative star-planct distance r, which, in 
general, is not constant due to the eccentricity e, is present. In 
this case, we have replaced it with the relative semi-major axis a 
because e = 0. 

^Note that eq. I|13|l automatically guarantees that the argument 
of the square root in oq. mill is positive. 



^ There is a wide variability in the possible values of fcj 
in hot Jupiters which can range from about 0.1 to 0.6 
l|Bodenheimer et al.ir200lf l. 

''Its precession is related to the change in tp; see, e.g., Iloridl|2007) l: 

Q 1 I2010I) . 
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Table 2 Relevant planet's parameters (jCoUier Cameron et al]|2010l ). We used mjup = 0.00189864 x 10^° kg, rjup = 
71.492 X 10^ m. The planet's angular speed ^p has been assumed equal to the K eplerian mean motion n = 27r/Pb, where 
the orbital period is Pb = 1.2198669 ± 0.0000012 d (jCoUier Cameron et aLlbOld). The values d isplayed for the rotational 



and tidal parameters gj? and qf are lower bounds corresponding to ( Collier Cameron et al.|[2010l ) 

plane t's mass nip. The value of the relative semi-major axis used in is a = 0.02555 ±0.00017 an ( Collier Cameron et al.l 



4.1 mjup for the 



I2OIOI ). Thefi gure for j'j has been obtained by using fcj — 0.3 for the planet's Love number. The value used for the flattening 
/p is the average of f^-^ and /p^x- 



Parameter (units) 


Numerical value 


rrip (mjup) 


< 4.1 


''p (^Jup) 


1.497 ±0.045 


?P (s-^) 


5.96147 X 10-^ 


q? 


0.00838747 


if 


-0.0251564 


3l 


2.09657 X 10^3 


fp 


0.00668828 


ftp 


0.277011 


Sp (kg m^ s^i) 


1.47243 X 10^9 



duration ( Miralda-Escudel 2002 ) . P reliminary numeri- 
cal an alyses have been performed by Ragozzine fc Woll 
( 2009[ ). 



with 



cos 



T < 



3.1 Classical and relativistic node precessions 

Concerning the tidal effects accounted for by qt, they 
do not secularly affect the node. Indeed, the tidal 
correction Utid to the pointlike two-bod y Newtonian 
gravitational potential o nly contains (|Sterndll939 : 
Ragozzine fc Wolj [20091 ) ■ so that when it is averaged 
over one orbital revolution the angle ip between the 
orbital angular momentum and the star/planet's spin 
axes, defined as '0 = when the two angular momenta 
are parallel and 1/' = tt when they are antiparallcl, does 
not appear in (J7tid); actually, the node secular pre- 
cession (} is proportional just to d (C/tid) /dip according 
to the Lagrange planetary equatio ns for the variation 
of th e Keplerian orbital elements (jMurrav fc Dermott 
2OOOI ). 

On the contrary, the centrifugal effects accounted for 
by Qr do affect the n odes because the potenti al [/rot 
contains the angle tp ( Ragozzine fc WoU 2009f ). The 
node precessions in the case of a two-body system with 
arbitrary masses niA and ms and rotational quadrupole 
mass moments and have been worked out by 



Barker fc O'Connelll (|l975l ) and IWexl (|l995il . In our 

case, the total precession is 



(22) 



(23) 



OP — 

qr ~ 2 



'rp -1 2 cos ipp 
V a I (l-e2)2 



(-5.880 X 10-")cosi/;p 

corresponds to 1.78 deg d^^, 
which is 9 X 10^ times larger than the node precession o f 
Mercury induced by the Sun's oblateness ( loriol 2005 ). 



Concerning the numerical value of O* quoted in eq. 



|), it has to be regarded as an upper limit computed 
for the largest possible value of cos^E**. Actually, the 
angle should not be confused w itljl A^, = 251.2 ±1.0 



deg ( Collier Cameron et al. 2010t ). which is, instead 



the projected angle between the star's spin axis and the 
orbital angular momentum defined in such a way that 
Aj, = deg if the sky-projected angular momenta are 



paral lel and A^ = tt if they are antiparallel (|Winn et al 
2OO9I) . In general, A^ can be derived from accurate 



radial- velocity data through the Rossiter-McLaughlin 



effect ([Anderson et al.ll2010[ ). but in the present case 



characterized by large stellar rotation its has been ob- 
tained from a line-profile analy sis in a far more direct 



and l ess model-dependent way (jCoUier Cameron et al, 
2OIOI) . The relation between ^f^,, J^,, A^,, i can be ob- 
tained as follows. From the spherical law of cosines 



A by ICoUier Cameron et al and 



^It is denoted a s 
IWinn et al.] (120091) . 
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Gellert et al. 


19891 


Zwillinger 


1995) 



witn 



cos A 



' cos C cos B + sin C sin B cos a, 



(24) 



the identifications A ^ 'i/^.^B tt — i,C — 7-/*,a— 7>A^, 
so that the segment a is due to the plane of the sky, the 
segment c is due to the orbital plan e and the segment b 
is the star's equatorial plane, yield ( Winn et al. l2009l ) 



cos = cos cos i + sin 7^, sin i cos A^,. 



(25) 



Given the determined values of i and A^,, eq. 
applied to the constraints on 7^, of cq. ([S|) , tells us that 
the minimum and the maximum values for ^'^^ are 



I* = ^min ^ *min = 91-30 deg, 

h = 97.20 dcg ^ ^-^^^^ = 108.93 deg, 



(26) 



i.e. the orbit of WASP-33b can be considered nearly 
polar. Consequently, we have 



|cos*^ 
I COS 5'* 



0.022, cos*^ < 0, 
0.324, cos*^ < 0. 



(27) 



The value h = h yields = 108.78 deg and cos = 
"0.32. Note that the constraints on of eq. (|26p are 
compatible with the observations. Indeed, since the 
time series of the residual average spectral line pro- 
file of HD 15082 during the transits unambiguously 
reveals that the motion of the planet is retrograde 
(jCollier Cameron et al.l 120101 ). it must be cosVP* < 0. 
Concerning i/jp, i.e. the angle between the planet's spin 
axis and the orbital angular momentum, it is not known 
and, at present, it seems that there is no way to measure 
ill. The largest possible value for corresponds to 
= Q deg. Indeed, it could be naively expected that 
the planet's rotation became synchronized and aligned 
with the orbital angular momentum on a rather short 
timescale; there might be a possibility that it gets into 
a Cassini state (jWinn fc Holmanll2005l) . 

The general relativistic gravitomagnetic (CM) spin- 
orbit effects in the case of a two-body system with 
arbitrary masses mA a nd ttt-b and spins 5 'a and 



5*6 have been derived by [Barker fc O'Connelll (|1975[ ). 
Damour fc Schafeil (|l988l ). IWexl (|l995l ). In moTcase. by 



posmg 



^ < 2.6 X 10-3, 



the total gravitomagnetic node precession is 



(28) 



(29) 



'(Collier Cameron A., private communication, June 2010). 



"GM 



"gm 



(1 + 3 ^ 2GS, 



(1 



(l_e2)3/2 



2Gap 
2a3(l-e2)3 



1 X 10-^^ s-i. 



(30) 



Contrary to the previously examined precessions due 
to the oblateness, the gravitomagnetic node precessions 
do not depend on the angles between the orbital angu- 
lar momentum and the bodies' spin axes. Note that 
1 X 10"^^ corresponds to 651 arcsec cty"^, which is 
3.25 X 10^ larger than the predicted solar Lense- Thirring 
precession of the node of Mercury ( loriol 2005 ). and 15 
times larger than the well-known gravitoelectric Ein- 
stein precession of the perihelion of Mercury in the field 
of the Sun. 

Table |3] summarizes the main precessional effects 
considered in this Section. The gravitomagnetic (star- 
induced) precession is 1 — 2 orders of magnitude smaller 
than both the oblateness precessions. The largest 
one is due to the star's oblateness; it overwhelms the 
planetary-induced precession by about 1 order of mag- 
nitude. Concerning the uncertainties in the nominal 
values of the dominant precessions fig^ , they can be con- 
servatively evaluated by linearly propagating the errors 
in the parameters entering their expressions. It turns 
out 



(5(f2*J < 0.106 
(5(f2Pj < 0.163 iiP^; 



(31) 



the uncertainties are of the same order of, or larger 
by 1 order of magnitude than the general relativis- 
tic gravitomagnetic precessions. Actually, they may 
be larger because we have kept the Love numbers ^2 
fixed to reference values. The previous considerations 
about the reciprocal orders of magnitude of the various 
node precessions hold for the upper bound of 17*^ due 
to I cos 4'*|niax, but it is easy to realize that, substan- 
tially, retain their validity also for | cos^^jmin yielding 
tl^^ = 5.169 X 10"" s-^ 

3.2 The effect of a third body 

The node precession induced by a distant third body X 
can be computed quite generally by using the Lagrange 



^"In the test particle limit for one of the two bodie s, we have the 
Lensc-Thirring precession ijLense fc ThirrinSllQlst ). Actually, in 
our case S^qj^^ is negligible, and the Lense-Thirring approximation 
is fully adequate for S^qj^. 
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Table 3 Node precessions due to the classical centrifugal oblateness and the general relativistic gravitomagnetic fields of 
both the star and the planet. The value of fij^ corresponds to | cos 3'*|max = 0.324. 



is-') 


is-') 


^GM ^) 


^GM is ^) 


< 7.6138 X 10-1" 


(-5.880 X 10-")cos?Ap 


1.53 X 10^12 


1 X 10-14 



planetary equations ( Murray &: DermottlboOO ) without 
making any a-priori assumptions on its location. The 
perturbing potential induced by X is ( Hogg et"al]ll991 ) 



Gmx ' 
-2^ 



3(r ■ h 



(32) 



where = fx/ fx is a unit vector pointing towards 
X. By denoting Ixjl-y, Iz the direction cosines of rjc , i.e. 
the components of Ix, we will express Ux as 



Ux = ^[r' ~3i^l. + yly + zh)'] 



(33) 



In working out the secular effects by X on the orbit of 
WASP-33b, eq. ([55)1 must be evaluated onto its unper- 
turbed Keplerian ellipse and averaged over one orbital 
revolution; in doing that, both rx and Ix will be rea- 
sonably assumed constant over one orbital period of 
WASP-33b, consistently with the hypothesis that if X 
exists, it must be quite distant. Useful formulas used 
in the calculation are 

X = r (cos fl cos u — cos sin sin u) , 
y = r (sin il cos u + cos cos sin u) , 
z = rsin'I'v.sinw, 

cos E~e 



cos/ = ^^1^, (34) 
sin / 
r 
dt 

where u == cj + / is the argument of latitude, / is the 
true anomaly, and E is the eccentric anomaly. The 
Gauss equatio ns for the v ariation of the Keplerian or- 
bital elements |RovII2005I ) straightforwardly yield 



1 — ecosi? ' 

a(l — ecosi?). 



<iO|(X)\ 
dt I / 



3Gmx CSC 'I't(l-e^ 
2rin 



,sin**- 



4- cos [ly cos Vl — Ix sin fi)] x 
X [Iz cos + sin X 



(35) 



An upper bound on the magnitude of the node preces- 
sion of WASP-33b due to a putative planet X can be 
naively obtained from 



\nx\ < 



for 



1 



3.18763 X 10^ 



mx = mjup. 



(36) 



(37) 



In Figure [T] we plot it as a function of its distance rx 
from 0.02 au to 10 au. It shows that the effect of the 
star's oblateness would be mimicked by a planet of one 
jovian mass at about rx ~ 0.1 — 0.15 au, while an ef- 
fect as large as the general relativistic gravitomagnetic 
precession would take place if X was at about 0.7 au. 
After 1 au, the impact of a jovian X would reduce to 
less than 10-^^ s-^, amounting to 10-^^ at rx ~ 10 



3.3 Measurability of the nodal precession 

Let us define as x reference direction the intersection 
between the plane of the sky, perpendicular to the line 
of sight, and the stellar equatorial plane, so that 



X = 



sin/* 



(38) 



in it p is the unit vector of the line of sight, oriented 
towards the observer. The line of nodes is the inter- 
section of the orbital plane with the star's equatorial 
plane, and forms an angle 51 relative to the reference x 
axis. The unit vector of the line of the node is 



S.xL 
sin** ' 



(39) 



where L is the unit vector directed along the orbital 
angular momentum. From eq. (|38p and eq. p9l) it 
follows 



sin vf* sin /* cos = (5* X L) ■ (5* X p). 



(40) 



Concerning the methods for measuring the node fJ, 
we mention that spectroscopic variations of the star 
during the transits can, in principle, be used to measure 



X {Ix sin — ly cos ^l)] 
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Maximum node precession due to a planet X Maximum node precession due to a planet X 




1 2 3 4 5 56789 10 

I'x (au) rx (au) 



Fig. 1 Approximate upper bound (3/2)(Gmx/ri)(l/n) of tlie node precession of WASP-33b due to a putative planet X 
for mx = mjup as a function of its distance for 0.02 au < rx < 10 au. 



the angle fl of the hne of nodeq^ ( Queloz et ah 2000 ). 
Recen tly, it has been pointed out by Fluri &: Berdyugina^ 

that the analysis of the polarization of the light 
scattered by the planets'atmospheres may allow to de- 
terminj^. among other things, as well within a 
few percent, according to preliminary numerical in- 
vestigations. Moreover, it turns out that thtj^ Q/F 
and U/ F polarized light curves are sensitive just to a 
time-variation of th e node , as shown by Figure 4 of 
Fluri fc Berdvugina ( 2010l ). Instead, the total degree 
of polarization P is not affected by dVl/dt. 

In the following, we will mainly focus our attention 
on the effects that a secular node variation has on the 
transit duration td for which relatively simple analytical 
expressions can be derived. 

3.4 The time variation of the transit duration 

The precession of the orbital plane changes the duration 
of the transit. For a circular orbit arbitrarily inclined 



Actually, the measurable quantity in iQueloz et al. I 1I2OOOI ') IS 
f! sin li, . 

With an ambiguity of 7r l|Fluri fc Berdvuginalboid ). 

^''They are the St okes coefficients defined by the Rayleigh scat- 
tering llFluri &: Berdvueina.i201Q) . 



to the line of sight the transit duration can be written 

as 



2{R, + rp) 
td = — cos 0, 



(41) 



where 5, the latitu de of the transit on the stellar di sk, 
is defined from0 (|Deedll998t lMiralda-Escudil2002[) 



sin (5 



a cos I 



(42) 



Intuitively, if i = 90 deg, the transit occurs along a 
diameter of the stellar disk, so that (5 = 0, as confirmed 
by eq. ([42|) : thus, eq. (|4T|) reduces to the intuitive 
fornJ^i 

2{R. + rp) 



td 



(43) 



^■*W hile iDeed lll998h -s ee hi s Figure 2-neglects the planet's ra- 
dius. |Miralda-E scud takes it into account. Note that 

iMiralda-Escud uses the angle between the orbital plane 

and the line of sight, and denotes it as a\ in t erms of our i, it 
is a = 7r/2 -|- i, so that sin a — > cosi. Moreover. iMiralda-Escudd 
adopts the letter 7 for the latitude of the transit 5. 

^^In this case, the transit duration is simply the time interval 
elapsed from the instant when the planet's disk touches the star's 
disk on one side of it to the instant when the planet's disk leaves 
the star's disk on the opposite side, thus traveling for a distance 
2R,i -I- 2rp at a speed given by its orbital velocity na. 
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Instead, for i = 0, no transit occurs; indeed, eq. (|42|) 
loses its meaning since it would imply sin (5 > 1. The 
quantity sin (5 can be thought as an adimensional im- 
pact parame ter, which is just one of the q uantities de- 
termined by Collier Cameron et al. ( 2010l) . By naming 
it &, they obtain 



0.155-0.218. 



(44) 



It should, now, be determined if the rate of change of 
td can be measured over a reasonable observing time 
baseline Atobs ~ 10 yr. Since the typical values of the 
precession periods is quite longer than 10 yr, only the 
time derivative dtd/dt is relevant. By differentiating eq. 
(|4ip with respect tc0 6 it is possible to obtain 



Tp) sin 5 / dS 



did 

dt na 
From eq. (gj) it follows 



dt 



dS_ 
'dt 



asm I 



di 



Tp) cos (5 \ dt 



which, substituted into eq. (|45p . yields 



dtd 
dt 



2 tan (5 sin i / di 
n \dt 



(45) 



(46) 



(47) 



From the spher ical law of cosines l|Gellert et al.l 11989 : 
ZwilhngeijflQoi ) 



cos B = sin C sin A cos b — cos C cos A 



(48) 



with the identificationg^l A — > ^'^,,5 tt - i,C ^ 
6 — > TT — 51, it turns out 



cos i = sin sin 7^, cos fl + cos cos 7^, . 



(49) 



From cq. (j49|) it is possible to obtain 51 as a function 
of A numerical investigation of it shows that f2 is 
close to 90 deg for the most likely values of By 
differentiating eq. with respect to fl one gets 



. . I di\ . ^ ■ T ■ f 
sm I \ -r ] = sm sm 1^ sm i I — — 
dt I \ dt 



(50) 



^^We do not take the derivative with respect to a because the 
secular variation of the semimajor axis vanishes for all the dy- 
namical effects considered. 

^'^In such a way, H results to be progradc with respect to the 
orbital motion, i.e. O follows it. 

^Miralda-Escudel \20Q± denotes /3 the angle between his mean 
plane-which in our case coincides with the star's equator-and the 
line of sigh t, so that /3 = 7r/2 — with this change, eq. (9) and 
eq. (12) of iMiralda-Escudj {iooj) agrees with our eq. l|49ll and 
eq. I|5 1t. It must also be noted that eq. (4) bv iMiralda-Escudel 
1120021 ) tells us that his ip and i coincide with our 



Finally, eq. (gZl) with eq. yields 

dtd 2 tan (5 sin sin sin 57 / dil 
dt n \ dt 



(51) 



To check intuitively the consistency of eq. (j49|) and of 
its consequences, let us consider a simplified situation 
in which the line of sight, the orbital angular momen- 
tum and the star's spin axis are mutually orthogonal. 
For example, we could imagine to see the orbital plane 
vertically in front of us, aligned with the stellar spin 
axis in an exactly polar configuration; actually, the real 
configuration of WASP-33b is not too far from the sce- 
nario described here. In this case, eq. (|49|) tells us that 
i = 51 = 90 deg, as expected. Moreover, according to 
our expectations, eq. (f50|) guarantees that the rate of 
change of i does not vanish, being just equal to that of 
57. The magnitude of the time variation of the transit 
duration induced by the star's oblateness through the 
node precession is shown in Figure [2] It is of the order 



Time variation of tlie duration transit due to stellar oblateness 



5 X lO-' 



4x10"' 



3 X lO-' 



axio-' 




Fig. 2 Time variation dtd/dt of tlie transit duration in- 
duced by the star's oblateness as a function of the angle 
7* between the line of sight and the stellar spin axis for 
11.22 deg < 7,^ < 168.77 deg. 



of 

dtd 
dt 



2-3 X 10" 



(52) 



over all the allowed range of values for 7,^. In Fig- 
ure [3] we plot (dtd/dt) cos ^pp due to the node preces- 
sion induced by the planet's centrifugal oblateness. It 
amounts to about —2 x 10~^ cos ipp for the most likely 
values of 7*. Figure |4] depicts the effect of the gravito- 
magnetic node precession on the transit duration. It is 
practically independent of 1^, , and amounts to 



dtd 
dt 



.4 X 10" 



(53) 



GM 



Thus, it is about 360 times smaller than the variation 
induced by the stellar oblateness and 24 times smaller 
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Time variation of the duration transit due to the planet's oblateness 



-5x 10- 




-2.5x10"' 



75 100 125 

I. (deg) 



Fig. 3 Time variation (dtd/dt) cos ipp of the transit du- 
ration induced by the planet's centrifugal oblateness as a 
function of the angle I, between the line of sight and the 
stellar spin axis for 11.22 deg < /* < 168.77 deg. 



Gravitomagnetic time variation of the duration transit 



8.41386 X Ur' 



8.41386 X Ur' 



75 100 125 

I. (deg) 



Fig. 4 Time variation dtd/dt of the transit duration in- 
duced by the gravitomagnetic node precession as a function 
of the angle 7, between the line of sight and the stellar spin 
axis for 11.22 deg < /* < 168.77 deg. 



than the maximum value of transit duration variation 
due to the planet's centrifugal oblateness. 

In order to evaluate the impact of a third body X on 
the variation of the transit duration, it is also necessary 
to derive eq. (|49)) with respect to ^E*^. Indeed, contrary 
to the centrifugal oblateness and general relativity, X 
causes a non- vanishing precession of the angle between 
the orbital and the star's equatorial plane as well. The 
perturbative Gauss equations yield 



dt I / 



^ (/. cos ^ + ly sin X 



\lz cos + sin X 



(Ix sin fl — ly cos fl)] 



Note that its magnitude is of the same order of mag- 
nitude than the node precession of eq. (|35|) . Thus, we 
have 



■ di 



(sin cos li, — cos 5*^, sin 7^, cos fl) 
sin^I/* sin/* sinil^, 



dt 



(55) 



so that 

dtd _ 
dt 



2 tan S 



[(sin^^cos/*- 



eos^I'j, sin/j, cos 17) 



dt 



(56) 



sinvl/*sin/*sinr2^] . 



It is useful to look at the multiplicative trigonometric 
factors of 'i>^, and (l in eq. A numerical investi- 

gation shows that the factor of Ct ranges from 1000 to 
4500 for the expected values of /* (Figure O, while the 
factor of ^fj, ranges from about —4000 to -1-4000 being 
zero close to /* = /* (Figure [6]) . 



Trigonometiicfactorof dil/dt in the transit duration variation 




(54) 



Fig. 5 Trigonometric factor of f2 in dtd/dt as a function 
of the angle J* between the line of sight and the stellar spin 
axis for 11.22 deg < /* < 168.77 deg. 

Concerning the accuracy with which a secular change 
of the transit duration may be measured over a typical 
time span At = 10 yr, it seems reasonable to a ssume 
a « 10"^ leveEl; see also iMiralda-Escudel |2002l ). who 
speculates about the possibility o f reaching a » 10~ ^ 
level, and the thorough analysis by Carter et al. ( 20081) . 
This would imply that the effect of the star's oblateness 
would be accurately measurable, while the detectabil- 
ity of the variation in the transit duration due to the 
planet's centrifugal oblateness strongly depends on the 
angle between the planet's spin axis and the orbital 



Collier Cameron A., private communication, June 2010. 
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Trigonometricfactorof d'P/dt in the transit duration variation 





-2000 
-4000 




25 50 75 100 125 150 

I. (deg) 



Fig. 6 Trigonometric factor of in dtd/dt as a function 
of the angle 7* between the line of sight and the stellar spin 
axis for 11.22 deg < J* < 168.77 deg. 

angular momentum. The general relativistic gravito- 
magnctic effect, instead, falls slightly below the mea- 
surability threshold. The same holds for a third body 
X with mx = JTijup located at about more than 0.7 au. 

4 Summary and conclusions 

The quadrupole mass moment and the proper angu- 
lar momentum 5*^, of the fast rotating, main sequence 
star WASP-33 are 1900 and 400 times, respectively, 
larger than those of the Sun. Thus, huge classical and 
relativistic non-Keplcrian orbital effects should affect 
the motion of the hot Jupiter harbored by WASP-33 
which has been recently discovered with photometric 
transit measurements. If measurable, they would yield 
important information on the physical properties of the 
star and the planet. 

Ill particular, the large inclination of the orbit 
of WASP-33b to the star's equator allows to consider 
the node precession VL. It could be determined from 
the polarimetry of the light curves, and the time vari- 
ation dtd/dt of the transit duration td- The node rate 
of WASP-33b induced by is 9 x 10^ times larger 
than the same effect for Mercury caused by the Sun's 
oblateness, while the general relativistic gravitomag- 
netic node precession-which is the only relativistic effect 
potentially measurable for WASP-33b if further analy- 
ses will confirm that its orbit is circular-is 3 x 10' times 
larger than the Lense-Thirring effect for Mercury due to 
the Sun's rotation. We also considered the effect of the 
centrifugal oblateness of the planet itself-thc tidal bulge 
raised on it by the star docs not affect the node-and of 
a putative distant third body X. A conservative eval- 
uation of the accuracy in measuring dtd/dt over 10 yr 
suggests a level w 10^*. The magnitudes of the induced 
time change in the transit duration are of the order of 



3 X X 10""^, 8 X 10"^ for J^, the planet's ro- 

tational oblateness and general relativity, respectively. 
They depend on the angle between the stellar spin 
axis and the line of sight which is, at present, unknown, 
and on the mass of the planet, for which only an up- 
per bound is nowadays available. A yet undiscovered 
planet X with the mass of Jupiter orbiting at at more 
than 1 au would induce a transit duration variation of 
less than 4 x 10"^. 

In conclusion, it appears that the star's oblateness 
can be measured at a percent accuracy from a 10-yr 
analysis of the time variations of the duration transit. 
The signals from the planetary rotational oblateness 
are smaller by 1 order of magnitude, at the least; ac- 
tually, they depend on the unknown angle '0p between 
the planet's equator and the orbital plane in such a way 
that they may turn out to be undetectable. The gen- 
eral relativistic gravitomagnetic signatures fall slightly 
below the measurability level. The presence in the sys- 
tem of a putative, distant third body, required in sev- 
eral theoretical schemes to explain the misalignment be- 
tween the stellar and orbital angular momenta, seems 
to be undetectable if it has the mass of Jupiter and 
orbits at more than 1 au from the star. 

The present analysis may be fruitfully repeated also 
if and when further planets moving along close inclined 
orbits around fast rotating stars will be discovered. 
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